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ABSTRACT 

Properties and structure of neutron stars are determined by the equation of 
state (EOS) of neutron-rich stellar matter. While the collective flow and particle 
production in relativistic heavy-ion collisions have constrained tightly the EOS 
of symmetric nuclear matter up to about five times the normal nuclear matter 
density, the more recent experimental data on isospin-diffusion and isoscaling in 
heavy-ion collisions at intermediate energies have constrained considerably the 
density dependence of the nuclear symmetry energy at subsaturation densities. 
Although there are still many uncertainties and challenges to pin down completely 
the EOS of neutron-rich nuclear matter, the heavy-ion reaction experiments in 
terrestrial laboratories have limited the EOS of neutron-rich nuclear matter in 
a range much narrower than that spanned by various EOSs currently used in 
astrophysical studies in the literature. These nuclear physics constraints could 
thus provide more reliable information about properties of neutron stars. Within 
well established formalisms using the nuclear constrained EOSs we study the mo- 
menta of inertia of neutron stars. We put the special emphasis on the component 
A of the extremely relativistic double neutron star system PSR J0737-3039. Its 
moment of inertia is found to be between 1.30 and 1.63 (xl0 45 g cm 2 ). Moreover, 
the transition density at the crust-core boundary is shown to be in the narrow 
range of p t = [0.091 - 0.093] (/m^ 3 ). 

Subject headings: dense matter — equation of state — stars: neutron — stars: 
rotation 



1. Introduction 



Neutron stars exhibit a large array of extreme characteristics. Their properties and 
structure are determined by the equation of state (EOS) of neutron-rich stellar matter at den- 
sities up to an order of magnitude higher than those found in ordinary nuclei, see, e.g. I Weber 
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( 119991 ) and lLattimer fc Prakashl (120041 ) . Therefore, the detailed knowledge about the EOS 
of neutron-rich nuclear matter over a wide range of densities is necessary for the study of 
neutron stars. For isospin asymmetric nuclear matter, various theoretical studies have shown 
that the energy per nucleon can be well approximated by 

E(p, 5) = E(p, 5 = 0) + E sym (p)5 2 + 0(5 4 ), (1) 

in terms of the baryon density p = p n + p p , the isospin asymmetry 5 = (p„ — p p ) / (p n + P P ) , the 
energy per nucleon in symmetric nuclear matter E(p, 5 = 0), and the bulk nuclear symmetry 
energy E sym (p). Here we report the results of a study on the moment of inertia and the core- 
crust transition density of neutron stars within well established formalisms in the literature 
using several EOSs constrained by the latest terrestrial heavy-ion reaction experiments. 

Presently, besides the possibilities of phase transitions into various non-nucleonic states 
the behavior of nuclear matter under extreme densities, pressures and/or isospin-asymmetry 
is still highly uncertain and relies upon, often, rather controversial theoretical predictions. 
This circumstance introduces corresponding uncertainties in the EOS of neutron-rich nuclear 
matter and thus limits our ability to understand many key issues in astrophysics. While as- 
trophysical observations can also limit the EOS of neutron-rich nuclear matter, terrestrial 
laboratories experiments provide complementary information and have their unique advan- 
tages. In this regard, it is especially interesting to mention that the collective flow and 
particle production in relativistic heavy-ion collisions have constrained the EOS of symmet- 
ric nuclear mat ter E(p, 5 = 0) up to about five tim es the normal nuclear matter density to a 



narrow range (IDanielewicz. Lacey. fc Lynch Il2002l ). However, there are still many challenges 



and uncertainties in pinning down precisely the EOS of neutron-rich nuclear matter. One of 



(Lattimer & Prakash 


2004; 


Steiner et al. 


2005; 


Chen et al. 


2007) 



E sym (p), see e.g. Refs. 
for recent reviews. To constrain the density dependence of the symmetry energy, many 
terrestrial nuclear experiments have been carried out recently or planned. Depending on 
the techniques used, some experiments are more useful for exploring the symmetry energy 
at low densities while others are more effective at high densities. For instance, heavy-ion 
reactions, especially those involving r adioactive beams, provide a unique means to probe the 



Esym i p) over a broad d ensity range (ILi et al.lll998t iLi fc Udo Schroederl l200ll ; iBaran et al. 



20051 ; IChen et al.l 120071 ). In fa ct, some significant progress has been made very recently by 



tudying the isos pin diffusion ( Shi fc Danielewic3 



Li fc Chen! 120051 ) and isoscaling (ITsang et al 



2003 



2001 



Tsang et al. 



2004 



Chen et al. 



Shetty. Yennello fc Souliotis 1 120071 ) in 



2005 



heavy-ion reactions at intermediate energies. The analysis of these phenomena based on 
transport theories of heavy-ion reactions and thermodynamical models of nuclear multifrag- 
mentation has limited the E sym (p) in a range much narrower than that spanned by various 
forms of the E sym (p) currently used in astrophysical studies in the literature. Moreover, 



- 3- 



the lower bound of the E sym (p) extracted from the heavy- ion reactions is consistent with 
the RMF prediction using the FSUGold interaction that can reproduce not only saturation 
proper ties of nuclear matt er but also structure properties and giant resonances of many finite 
nuclei (IPiekraewicz 1 120071 ). 



It is also well known that the sizes of neutron skins in heavy nuclei are sensitive to the 
symmetry energy at subsaturation densities, see, e.g., Refs. ( Brown II2OOO i Horowitz fc Piekarewicz) 
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Steiner et al.112005 



Todd-Rutel fc Piekarewicz 



ness obtained using hadronic probes are not accurate enough to constrain significantly the 
symmetry energy. Interestingly, the parity radius experiment (PREX) at the Jefferson Labo- 
ratory aiming to measure the neutron radius in 208 Pb via parity viola ting; electron scattering 
(Jefferson Laboratory Experiment E-00-003) (IHorowitz et al.l l200ll ) hopefully will provide 
much more precise data and thus constrain the symmetry energy at low densities more 
tightly in the near future. On the other hand, at supranor mal densities, a nu mber of po 



tential probes of the symmetry energy have been proposed (ILi et al.l 119971 ; iLil 12000 . 12002 



Li et al.l 


1998: 


Chen et al. 


2007) 



behavior of the symmetry energy with high energy radioactive beams have been planned at 
the CSR/Lanzhou, FAIR/GSI, RIKEN and the NSCL/MSU. 

While the EOS of neutron-rich nuclear matter has not been completely determined yet, 
it is still very interesting to examine astrophysical implications of the EOS constrained by 
the latest terrestrial laboratory experiments mentioned above. Global properties of spheri- 
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i|). However, properties of (rapidly) rotating neutron stars have 



const 


ructed only by several research groups with various degree of approximation ('. 


lartle 


1967; 


Hartle & Thorne 


19681: Friedman et al. 1986: 


Bombaci et al 12000 : Lattimer et al. 


1990; 


Komatsu et al. 1989: Cook et al. 1994: 


Stergioulas k, Friedman 


19951. 19981: Bonazzola et al. 


1993. 


1998; 


Weberl 199J 


); i 


Vnsorg et al. 


2002) (see 


Stergioulas ( 


2003) for a review). In 


1 a re- 
tasses, 


cent work ( 


Krastev et al. 


2008) we have reported predictions on the gravitational n 



radii, maximal rotational (Kepler) frequencies, and thermal properties of (rapidly) rotating 
neutron stars. In this work, using the nuclear constrained EOSs we calculate the momenta 
of inertia for both spherically-symmetric (static) and (rapidly) rotating neutron stars us- 
ing well established formalisms in the literature. Such studies are important and timely as 
they are related to the astrophysical observations in the near future. In particular, the mo- 
ment of inertia of pulsar A in the extremely relativistic neutron star binary PSR J0737-3039 



-4- 



( IBurgay et al.l 120031 ) m ay be determined i n a few years through detailed measurements of 
the periastron advance (IBejger et al.ll2005l ). 



2. The equation of state of neutron-rich nuclear matter constrained by recent 

data from terrestrial heavy-ion reactions 

In this section, we first outline the theoretical tools one uses to extract information about 
the EOS of neutron-rich nuclear matter from heavy-ion collisions. We put the special empha- 
sis on exploring the density- dependence of the symmetry energy as the study on the EOS of 
symmetric nuclear matter with heavy-ion reactions is better k nown to the astrophysical com- 



muni t y and it has been e xtensively reviewed, see e. g., Refs. ( jDanielewicz. Lacey. fc Lynch 



20021 ; ISteiner et al.l 120051 ; lLattimer fc Prakashl 120071 ) for recent reviews. We will then sum- 
marize the latest constraints on the density dependence of the symmetry energy extracted 
from studying isospin diffusion and isoscaling in heavy-ion reactions at intermediate ener- 
gies. Finally, we address the question of what kind of isospin-asymmetry, especially for dense 
matter, can be reached in heavy-ion reactions. 

Heavy-ion reactions are a unique means to create in terrestrial laboratories dense nu- 
clear matter similar to those found in the core of neutron stars. Depending on the beam 
energy, impact parameter and the reaction system, various hadrons and/or partons may 
be created during the reaction. To extract information about the EOS of dense matter 
from heavy-ion reactions requires careful modelling of the reaction dynamics and selec- 
tion of sensitive observables. Among the available tools, computer simulations based on 
the B oltzmann-Uehling-Uhlenbeck ( buu) transport theory have been ve ry useful, see, e.g., 
Refs. (IBertsch fc Das Gupta 1 1 19881 ; iDanielewicz. Lacey. fc Lynch 1 12002| ) for reviews. The 
evolution of the phase space distribution function fi(r,p,t) of nucleon i is governed by both 
the mean field potential U and the collision integral Illusion via the BUU equation 



dt 



+ V P U ■ V r / 8 - V r U ■ V p fi = I, 



collision- 



(2) 



Normally, effects of the collision integral I C oiUsion y i a both elastic and inelastic channels in- 
cluding particle productions, such as pions, are modelled via Monte Carlo sampling using 
either free-space experimental d ata or calculated in-medium cross sections for the elemen- 
tary hadron-hadron scatterings (IBertsch fc Das Gupta Ill988l ). Information about the EOS 
is obtained from the underlying mean-field potential U which is an input to the transport 
model. By comparing experimental data on some carefully selected observables with trans- 
port model predictions using different mean-field potentials corresponding to various EOSs, 
one can then constrain the corresponding EOS. The specific constrains on the density de- 
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pendence of the nuclear symmetry ener gy that we are us ing in this work were obtained 
by analyzing the isospin diffusion data (ITsang et al.l I2004T) within the IBUU04 version of 



(Lietal 



Das et al. 



2004). In this model, an 



20031 ) is used. With this 



an isospin and momentum dependent transport mode 
isospin and momentum-dependent interaction (MDI) ( 
interaction, the potential energy density V(p, T, 5) at total density p, temperature T and 
isospin asymmetry S is 



V(p,T,S) 



CT+l 



AupnPp Ai 2 2 . B p' 



Po 

Cry 



xS 2 ) 



Po 




ap(1P l + (p-^)2/ A 2 



(3) 



In the mean field approximation, Eq. ([3]) leads to the following single particle potential for 
a nucleon with momentum p and isospin r 



U T (p,T,6,p,x) = A^x^ + Mx^ + bI^] (1 

Po Po \PoJ 



x6 2 



ITX- 



B p 



a-l 



a 



7- SP- 

1 Po 



t=T, — T 



Po 



ft{r,p 



1 + (p-^) 2 /A 2 ' 



(4) 



where r = 1/2 (—1/2) for n eutrons (protons '), x, A u (x), Ai(x), B, C T)T ,C T - T , a, and A are 
all parameters given in Ref. (IDas et al.ll2003l ). The last two terms in Eq. (J3J) contain the mo- 
mentum dependence of the single-particle potential, including that of the symmetry poten- 
tial if one allows for different interaction strength parameters C Tj -_ T and C Tjr for a nucleon of 
isospin r interacting, respectively, with unlike and like nucleons in the background fields. It is 
worth mentioning that the nucleon isoscalar potential estimated from U isosca i ar ~ (U n + U p )/2 
agre es with the pre diction of variational many-body calculations for symmetric nuclear mat- 
ter (I Wiringal Il988l ) in a broad density and momentum range (ILi et al.l 12004 ). Moreover, 
the EOS of symmetric nuclear matter for this interaction is consistent with that extracted 
from the available data on collective flow and particle production in relativistic heavy-ion 



collisions up to five times the normal nuclear matter (IDanielewicz. Lacey. fc Lynch 1 12002 



Krastev et al. 1 120081 ). On the other hand, the corresponding isovector (symmetry) potential 
can be estimated from U sym ~ (U n — U P )/2S. At normal nuclear matter density, the MDI 
symmetry potential agrees very well with the Lane potential extracted from nucleon-nucleus 
and ( n,p) charge ex change reactions available for nucleon kinetic energies up to about 100 
MeV ( ILi et al.l 120041 ) . At abnormal densities and higher nucleon energies, however, there is 
no experimental constrain on the symmetry potential available at present. 



The different x values in the MDI interaction are introduced to vary the density depen- 
dence of the nuclear symmetry energy while keeping other properties of the nuclear equation 
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Fig. 1. — The density dependence of the nuclear s ymmetry energ y for different values of the 
parameter x in the MDI interaction. Taken from (ILi et al.ll2005l ). 



of state fixed. Specifically, choosing the incompressibility K$ of cold symmetric nuclear mat- 
ter at saturation density p to be 211 MeV leads to the dependence of the parameters A u 
and Ai on the x parameter according to 



AJx) = -95.98 -x- 



2B 

'<7+T 



Mx) 



-120.57 + a; 



IB 
a+1' 



(5) 



with B = 106.35 MeV. 



With the potential contribution in Eq. [3] and the well-known contribution from nucleon 
kinetic energies in the Fermi gas model, the EOS and the symmetry energy at zero tempera- 
ture can be easily obtained. As shown in Fig. (TJ adjusting the parameter x leads to a broad 
range of the density dependence of the nuclear symmetry energy, similar to those predicted 
by various micros copic and /or phen o meno logical many-body theories. As demonstrated by 
Li &: Chen! (120051 ) and ILi fc Steinerl (120061 ). only equations of state with x between -1 and 
have symmetry energies in the sub-saturation density region consistent with the isospin 
diffusion data and the available measurements of the skin thickness of 208 Pb using hadronic 
probes. Moreover, it is interesting to note that the symmetry energy extracted very recently 
from the iso scalinK analyses of heavy- ion react ions is consistent with the MDI calculation 
using x = (jShetty. Yennello fc Souliotis 1 120071 ). The E x ym (p) with x = is also consistent 



with the RMF prediction using the FSUGold interaction (IPiekraewicz 1 120071 ). We thus con- 
sider only the two limiting cases with x = and x = — 1 as boundaries of the symmetry 
energy consistent with the available terrestrial nuclear laboratory data. 
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Table 1: The parameter s F (MeV) G , K sym (MeV), L (MeV), and fC OS2/ (MeV) for different 
values of x. Taken from ()Chen et al.ll2005l ). 



X 



G 



K 



sym 



K, 



a. si/ 



1 107.232 1.246 -270.4 

129.981 1.059 -88.6 

-1 3.673 1.569 94.1 

-2 -38.395 1.416 276.3 



16.4 -368.8 

62.1 -461.2 

107.4 -550.3 

153.0 -641.7 



To ease comparisons with other models in the literature, it is useful to parameterize 
the E sym (p) from the MDI interaction and list its characteristics. Within phenomenological 
models it is customary to sep arate the symmetry energy into the kinetic and potential parts, 
see, e.g. ( jPrakash et al.lll9881 ). 



E sym {p) = (2 2 / 3 



,2/3 



(6) 



With the MDI interaction, the potential part of the nuclear symmetry energy can be well 
parameterized by 

E^ m (p) = F(x)p/ Po + (18.6 - F(x))(p/p ) GW , (7) 
with F(x) and G(x) given in Table [1] for x — 1, 0, —1 and —2. The MDI parameteriza- 



tions for the E^J^p) is similar but significantly different from those used by iPrakash et al. 
(119881 ). Also shown in Table [T] are other characteristics of the symmetry energy, including 
its slope parameter L and curvature parameter K sym at p , as well as the isospin-dependen t 
part i\asy of the isobaric incompressibility of asymmetric nuclear matter (IChen et al.l 120051 ). 



The symmetry energy in the subsaturation density region with x=0 and -1 can be roughly 
approximately by E sym (p) ps 31.6(p/p ) ' 69 and E sym (p) ps 31.6(p/p ) 105 , respectively. 

The MDI EOS has been recently applied to constrain the mass-radius correlation s of 
both static and rapidly rotating neutron stars (ILi &: Steinerll2006l ; iKrastev et al. 1 120081 ). In 



addition, it has been also used to constrain a possible time variation of the gravitational 
constant G (IKrastev fc Lil 120071 ) via the gravitochemical heating formalism developed by 
Jofre et al.l (120061 ). For comparisons, in this work we apply also EOSs fr om variational cal- 
culations with the Al8+Sv+UIX* interaction fAPRdAkmal et al.l (119981). and recent Dirac- 



Brue ckner-Hartree-Fock (DBHF) calculations (lAlonso fc Sammarruca 



2003; Krastev fc Sammarruca 



2006[) (DBHF+Bonn B) with Bonn B One-Boson-Exchange (OBE) potential JMachleidt 



19891 ). Below the baryon density of approximately 0.07/m 3 the equations of state are sup- 
plemented by a crustal EOS, which is more suitable for the low density regime. Namely, 
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Table 2: Saturation proper ties of the nuclear EO Ss (for symmetric nuclear matter) employed 
in this work. Taken from (IKrastev et al. II2008I ). 



EOS 


po(fm- 3 ) 


E s (MeV) 


K(MeV) 


e sym (p )(MeV) 


m*(po)/m 


MDI(x=0) 


0.160 


-16.08 


211.00 


31.62 


0.67 


MDI(x=-l) 


0.160 


-16.08 


211.00 


31.62 


0.67 


APR 


0.160 


-16.00 


266.00 


32.60 


0.70 


DBHF+Bonn B 


0.185 


-16.14 


259.04 


33.71 


0.65 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities at the nuclear saturation density: saturation (baryon) density; energy-per-particle; com- 
pression modulus; symmetry energy; nucleon effective mass to average nucleon mass ratio (with 
m = 938.926MeF c" 2 ). 



we apply the EOS by 



Pethick. Ravenhall &: Lorenzl (119951 ) for the inner crust and the one 



by lHaensel fc Pichonl (119941 ) for the outer crust. At the high est densities we assume a con 



tinuous functional for the EOSs employed in this work. (See (IKrastev fc Sammarrucall2006l ) 
for a detailed description of the extrapolation procedure for the DBHF+Bonn B EOS.) The 
saturation properties of the nuclear equations of state used in this paper are summarized in 
Table [2 

What is the maximum isospin asymmetry reached, especially in the supra-normal den- 
sity regions, in typical heavy-ion reactions? How does it depend on the symmetry energy? 
Do both the density and isospin asymmetry reached have to be high simultaneously in or- 
der to probe the symmetry energy at supra-normal densities with heavy-ion reactions? The 
answers to these questions are important for us to better understand the advantages and 
limitations of using heavy-ion reactions to probe the EOS of neutron-rich nuclear matter 
and properly evaluate their impacts on astrophysics. 

To answer these questions we first show in Fig. [2] the central baryon density (upper 
window) and the average {n/p) p > po ratio (lower window) of all regions with baryon densi- 
ties higher than po i n the reaction of 132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon 
and an impact parameter of 1 fm. It is seen that the maximum baryon density is about 
2 times normal nuclear matter density. Moreover, the compression is rather insensitive to 
the symmetry energy because the latter is relatively small compared to the EOS of sym- 
metric nuclear matter around this density. The high density phase lasts for about 15 fm/c 
from 5 to 20 fm/c for this reaction. It is interesting to see in the lower window that the 
isospin asymmetry of the high density region is quite sensitive to the density dependence 
of the symmetry energy used in the calculation. The soft (e.g., x — 1) symmetry energy 
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Fig. 2. — Central baryon density (upper panel) and isospin asymmetry (lower panel) of high 
density region in the reaction of 132 Sn + 124 Sn at a beam ene rgy of 400 MeV/nucleon and an 
impact parameter of 1 fm. Taken from Ref. (ILi et al.ll2005l ). 



leads to a significantly higher value of {n/p) p > po than the stiff one (e.g., x = —2). This 
is consi stent with the well-known isospin frac tionation phenomenon in asymmetric nuclear 
matter jMuller fc Serot Hl99fitEi fc Ko Ill997h . Because of the E sym (p)5 2 term in the EOS of 
asymmetric nuclear matter, it is energetically more favorable to have a higher isospin asym- 
metry 8 in the high density region with a softer symmetry energy functional E sym (p). In the 
supra-normal density region, as shown in Fig. [TJ the symmetry energy changes from being 
soft to stiff when the parameter x varies from 1 to -2. Thus the value of {n/p) p > Po becomes 
lower as the parameter x changes from 1 to -2. It is worth mentioning that the initial value 
of the quantity (n/p) p > po is about 1.4 which is less than the average n/p ratio of 1.56 of 
the reaction system. This is because of the neutron-skins of the colliding nuclei, especially 
that of the projectile 132 Sn. In the neutron-rich nuclei, the n/p ratio on the low-density 
surface is much higher than that in their interior. Also because of the E sym (p)5 2 term in the 
EOS, the isospin-asymmetry in the low density region is much lower than the supra- normal 
densi ty region as long a s the symmetry increases with density. In fact, as shown in Fig. 2 of 
Ref. ( Li fc Steinerl 120061 ). the isospin-asymmetry of the low density region can become much 
higher than the isospin asymmetry of the reaction system. 
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It is clearly seen that the dense region can become either neutron-richer or neutron- 
poorer with respect to the initial state depending on the symmetry energy functional E sym (p) 
used. As long as the symmetry energy increases with the increasing density, the isospin 
asymmetry of the supra-normal density region is always lower than the isospin asymmetry 
of the reaction system. Thus, even with radioactive beams, the supra-normal density region 
can not be both dense and neutron-rich simultaneously, unlike the situation in the core of 
neutron stars, unless the symmetry energy starts decreasing at high densities. The high 
density behavior of the symme try energy is probably among the most u ncertain properties 
of dense matter as stressed by (IKutschera et al. Ill993l ; iKutschera Il2000l ). Indeed, some pre- 
dictions show that the symmetry energy can decrease with increasing density above certain 
density and may even finally becomes negative. This extre me behavior was first predicted 



by some micro s copic many-body theories, see e.g., Refs. (IPandharipande &: Garde I Il972 
Wiringal Il988al ; iKrastev &: Sammarrucal 120061 ) . It has also been shown that the symme- 



try energy can become negative at various high densi ties within the Hartree-Fock approach 
using the original Gogny force (jChabanat et al.lll997l ). the density-dependent M3Y interac- 
tion (IKhoa et al.lll996t iBasu et al.l 120061 ) and about 2/3 of the 87 Skyrme interactions that 
have been widely used in the literature (IStone et al. 1120031 ). The mechanism and physical 
meaning of a negative symmetry energy are still under debate and certainly deserve more 
studies. 

Isospin effects in heavy- ion reactions are determined mainly by the E sym (p)S 2 term in 
the EOS. One expects a larger effect if the isospin-asymmetry is higher. Thus, ideally, one 
would like to have situations where both the density and isospin asymmetry are sufficiently 
high simultaneously as in the cores of neutron stars in order to observe the strongest effects 
due to the symmetry energy at supra-normal densities. However, since it is the product of the 
symmetry energy and the isospin-asymmetry that matters, one can still probe the symmetry 
energy at high densities where the isospin asymmetry is generally low with symmetry energy 
functionals that increase with density. Therefore, even if the high density region may not 
be as neutron-rich as in neutron stars, heavy-ion collisions can still be used to probe the 
symmetry energy at high densities useful for studying properties of neutron stars. 
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The moment of inertia of neutron stars 



Employing the EOSs described briefly in Section 2, we compute the neutron star moment 
of inertia w ith the RN^_ code developed and made available to the public by Nikolaos 
Stergioulas (jStergioulas fc Friedman! Il995l ). The code solves the hydrostatic and Einstein's 
field equations for mass distributions rotating rigidly under the assumption of stationary and 
axial symmetry about the rotational axis, and reflec tional symmetry about the equatorial 
plane. RNS calculates the angular momentum J as (jStergioulasI 120031 ) 



where T^ u is the energy-momentum tensor of stellar matter 

T" v = (e + P)itV + Pg^, (9) 

£^ is the Killing vector in azimuthal direction reflecting axial symmetry, and dV = y/—gd 3 x 
is a proper 3- volume element (g = det(g a p) is the determinant of the 3-metric). In Eq. ([9]) 
P is the pressure, e is the mass-energy density, and is the unit time-like four-velocity 
satisfying u^u^ = — 1. For axial-symmetric stars it takes the form = tt*(l, 0, 0, Q), where 
Q is the star's angular velocity. Under this condition Eq. (jSj) reduces to 



J = j(e + P)u\g H u <t> + g^ t u t )^gd 3 x 



(10) 



It should be noted that the moment of inertia cannot be calcula ted directly as an integral 
quantity over the source of gravitational field (jStergioulasI 120031 ) . In addition, there exists 
no unique generalization of the Newtonian definition of the moment of inertia in General 
Relativity and therefore I = J/Q is a natural choice for calculating this important quantity. 

For rotational frequencies much lower than the Kepler frequency (the highest possible 
rotational rate supported by a given EOS), i.e. vjv\. « 1 (u = ft/ (2n)), the deviations from 
spherical symmetry are very small, so that the moment of inertia can be approximated from 
spherical stellar mode ls. In what follows we review briefly this slow-rotation approximation, 
see e.g. ( iHartld 119671 ). In the slow-rotational limit the metric can be written in spherical 
coordinates as (in geometrized units G = c = 1) 



ds 1 



- e ^(r) dt 2 + 1 



2m(r) 



dr 2 - 2ur 2 sin 2 9dtd(j) + r 2 (d6 2 + sin 2 



1 Thanks to Nikolaos Stergioulas the RNS code is available as a public domain program at 



http : / / www. gravity, phys . uwm.edu/rns / 



In the above equation m(r) is the total gravitational mass within radius r satisfying the 
usual equation 

dm(r) 2 



dr 



4ire(r)r' 



(12) 



and u)(r) = (d<p/dt)zAMO is the Lense-Thirring angular velocity of a zero-angular-momentum 
observer (ZAMO). Up t o first order in u all m etric functions remain spherically symmetric 
and depend only on r (IMorrison et al.l l2004t ) . In the stellar interior the Einstein's field 
equations reduce to 



d(p{r) 
dr 



m(r) 



1 + 



AiTr 3 P(r) 



m r 



1 - 



2m(r) 



~i -i 



(r <R 



star I 



and 



1 d 

r 3 dr 



4 v ,duj(r) 
r jir)—— 
dr 



dr 



(r <R 



star ) i 



(13) 



(14) 



with uj = Q — u the dragging angular velocity (the angular velocity of the star relative to a 
local inertial frame rotating at u) and 



J 



1 - 



2m(r] 



1/2 



Outside the star the metric functions become 

2M 



D 2</> 



(r > R star ) 



and 



2J 



(r > R 



star ) i 



(15) 



(16) 



(17) 



where M = m(r = R) = Arc e(r')r' 2 dr' is the total gravitational mass and R is the stellar 
radius defined as the radius at which the pressure drops to zero (P(r = R) = 0). At the 
star's surface the interior and exterior solutions are matched by satisfying the appropriate 
boundary conditions 

«-f(fL (18) 

and 



bn 
2 



2M 



(19) 



The moment of inertia I = J/Q then can be computed from Eq. (fTOj) . With Q = u^/ii* and 



retaining only first order terms in uj and Q, the moment of inertia reads (Morrison et al. 



2004 : lLattimer & Prakashll2000f ) 



13 



This slow-rotation approximation for the neutron-star moment of iner tia neglects deviation s 
from spherical symmetry and is independent of th e angular velocity Q (IMorrison et al.ll2004j ). 
For neutron stars with masses greater than 1M lLattimer fc Schutzl (120051 ) found that, for 
slow-rotations, the momenta of inertia computed through the above formalism (Eq. (|20p ) 
can be approximated very well by the following empirical relation: 



(0.237 ±0.008)Mi? 2 



, Mkm ( Mkm\ 4 

1 + 42 m^ + 9( WJ 



(21) 



The above equation is shown (lLattimer fc Schuta 120051 ) to hold for a wide class of EOSs 
except for ones with appreciable degree of softening, usually indicated by achieving a maxi- 
mum mass of ~ 1.6M or less. Since none of the EOSs employed in this paper exhibit such 
pronounced softening, Eq. ( 12T|) is a good approximation for the momenta of inertia of slowly 
rotating stars. 




Fig. 3. — (Color online) Mass-radius relation for static and maximally rotating neutron stars. 



ines correspond to st atic and broken lines to maximally rotating stellar models. Taken 



Solid 



form (IKrastev et al. 1120081 ). 



We calculate the neutron star moment of inertia applying several nucleonic EOSs (see 
Section 2) considering both slow- and rapid-rotation regimes. In Fig. [3] we show stellar se- 
quences computed with the RNS code for spherical and maximally rotating models. As seen 
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in the figure, rapid rotation alters significantly the mass-radius relation of rapidly rotating 
stars (with respect to static configurations). Generally, for a given EOS it increases the max- 
imum possible mass by ~ 15%, while reducing/increasing the polar/circumferential radius 
by several kilometers, leading to an overall oblate shape of the rotating star. The degree to 
which the neutron star properties are impacted by rapid rotation depends on the details of 
the EOS: it is greater for models fro m stiffer EOS which pr oduce less centrally condensed and 
gravitational bound neutron stars ([Friedman et al.lll984l ). In view of these considerations, 
one should expect similar changes in the moment of inertia of rapidly rotating neutron stars 
(with respect to static models). We address these and other implications next. 



4.1. Slow rotation 



If the rotational frequency is much smaller than the Kepler frequency, the deviations 
from spherical symmetry are negligible and the moment of inertia can be c alculated applying 
the slo w-rotation approximation discussed briefly in Section 3. For this case lLattimer fc Schutz 
(120051 ) showed that the moment of inertia can be very well approximated by Eq. (l2Tj) . In 
Fig. H] we display the moment of inertia as a function of stellar mass for slowly rotating neu- 
tron stars as computed with the empirical relation ( !2~lT) . As shown in Fig. [3j above ~ 1.0M Q 
the neutron star radius remains approximately constant before reaching the maximum mass 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 

M(M @ ) 



Fig. 4. — (Color online) Total moment of inertia of neutron stars estimated with Eq. (l2Tl) . 
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x=-1 

DBHF+Bonn B 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 

M(M @ ) 



Fig. 5. — (Color online) The moment of inertia scaled by M 3//2 as a function of the stellar 
mass M. The shaded band illustrates a 10% error of hypothetical I/M 3 ^ 2 measurement of 
50 km 2 M~ 1/2 . Th e erro r bar shows the specific case in which the mass is 1.34M (after 



flLattimer fc Schutzl 120051 )1 



supported by a given EOS. The moment of inertia (/ ~ MR 2 ) thus increases almost linearly 
with stellar mass for all models. Right before the maximum mass is achieved, the neutron 
star radius starts to decrease (Fig. [3]), which causes the sharp drop in the moment of inertia 
observed in Fig. @J Since I is proportional to the mass and the square of the radius, it is 
more sensitive to the density dependence of the nuclear symmetry energy, which determines 
the neutron star radius. Here we recall that the x = — 1 EOS has much stiffer symmetry 
energy (with respect to the one of the x = EOS), which results in neutron star models 
with larger radii and, in turn, momenta of inertia. For instance, for a "canonical" neutron 
star (M = 1.4M ), the difference in the moment of inertia is more than 30% with the x = 
and the x = —1 EOSs. In Fig. [5] we take another view of the moment of inertia where I is 
scaled by M 3//2 as a function of the stellar mass (after ( Lattimer &: Schutzl 120051 )). 



The discovery of the extremely relativistic binary pulsar PSR J0737-3039 A,B provides 



an un p recedented opportunity to test General Relativity and physics of pulsars (jBurgay et al. 



20031 ). lLattimer fc Schutzl (120051 ) estimated that the moment of inertia of the A component 
of the system should be measurable with an accuracy of about 10%. Given that the masses of 
both stars are already accurately determined by observations, a measurement of the moment 
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Table 3: Numerical results for PSR J0737-3039A (M A = 1.338M©, v A = 44.05Hz). 



EOS 




e c (x W 14 g cm 3 ) 


R eq (km) 


J(xl0 45 # cm 2 ) 


I LS {xl0 45 g cm 2 ) 


MDI(x= 


=-1) 


7.04 


13.75 (13.64) 


1.63 


1.67 


DBHFH 


-Bonn B 


7.34 


12.56 (12.47) 


1.57 


1.43 


MDI(x= 


=0) 


9.85 


12.00 (11.90) 


1.30 


1.34 


APR 




9.58 


11.60 (11.52) 


1.25 


1.26 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities: central mass-energy density, equatorial radius (the numbers in the parenthesis are the 
radii of the spherical models; the deviations from sphericity due to rotation are ~ 1%), total 
moment of inertia, total moment of inertia I LS as computed with Eq. (|2ip . 



of inerti a of even one neutron st ar could have enormous importance for the neutron star 
physics ( ILattimer &: Schutzl 120051). (The signifi cance of such a measurement is illustrated in 
Fig.[5J As pointed by ILattimer &: Schuta (120051 ). it is clear that very few EOSs would survive 
these constraints.) Thus, theoretical predictions of the moment of inertia are very timely. 
Calculation s of the moment of iner tia o f pulsar A {Ma = 1. 338M , v& = 4:4.05Hz) have been 
reported by iMorrison et al.l (120041 ) and lBejger et al.1 (120051 ). In Table [3] we show the moment 
of inertia and (other selected quantities) of PSR J0737-3039A computed with the RNS code 
using the EOSs employe d in this study. Our results with the APR EOS are in very good 
agreem ent with those by IMorrison et al.1 (120041 ) (/ APi? = 1.24 x 10 45 g cm 2 ) and lBejger et al. 
( 120051 ) {I APR = 1.23 x 10 45 g cm 2 ). In the last column of Table [3] we also include results 
computed with the empirical relation (Eq. (12T1) ). From a comparison with the results from 
the exact numerical calculation we conclude that Eq. (|2l|) is an excellent approximation for 
the moment of inertia of slowly-rotating neutron stars. (The average uncertainty of Eq. (12T1) 
is ~ 2%, except for the DBHF+BonnB EOS for which it is ~ 8%.) Our results (with the 
MDI EOS) allowed us to constrain the moment of inertia of pulsar A to be in the range 
I = (1.30 - 1.63) x 10 45 (g cm 2 ). 



4.2. Rapid rotation 



In this subsection we turn our attention to the moment of inertia of rapidly rotating neu- 
tron stars. In Fig. [6] we show the moment of inertia as a function of stellar mass for neutron 
star models spinning at the mass-shedding (Kepler) frequency. The numerical calculation is 
performed with the RNS code. We observe that the momenta of inertia of rapidly rotating 
neutron stars are significantly larger than those of slowly rotating models (for a fixed mass). 
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Fig. 7. — (Color online) Total moment of inertia as a function of stellar mass for models 
rotating at 716Hz (upper frame) and 1122 Hz (lower frame). 
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Fig. 8. — (Color online) Total moment of inertia as a function of rotational frequency for 



stellar models with mass M = I AM, 



©■ 



This is easily understood in terms of the increased (equatorial) radius (Fig. [3]). 



We also comp ute the momenta o f inertia of neutron stars rotating at 716 (IHessels et al. 



20061 ) and 1 122Hz (IKaaret et al.ll2007l ) which are the rotational frequencies of th e fastest pul 



sars o f today. The numerica l results are presented in Fig. [TJ As demonstrated by lBejger et al. 



( 120071 ) and most recently by lKrastev et al. I (120081 ). the range of the allowed masses supported 
by a given EOS for rapidly rotating neutron stars becomes narrower than the one for static 
configurations. The effect becomes stronger with increasing frequency and depends upon 
the EOS. This is also illustrated in Fig. [7J particularly in the lower panel. Additionally, the 
moment of inertia shows increase with rotational frequency at a rate dependent upon the 
details of the EOS. This is best seen in Fig. [S] where we display the moment of inertia as a 
function of the rotational frequency for stellar models with a fixed mass (M = 1.4M ). The 
neutron star sequences shown in Fig. [8] are terminated at the mass-shedding frequency. At 
the lowest frequencies the moment of inertia remains roughly constant for all EOSs (which 
justifies the application of the slow- rotation approximation and Eq. (12TT) ). As the stellar 
models approach the Kepler frequency, the moment of inertia exhibits a sharp rise . This 
is attributed to the large increase of t he circumferent i al rad ius as the star approaches the 
"mass-shedding point". As pointed by lFriedman et al.l (119841 ). properties of rapidly rotating 
neutron stars display greater deviations from those of spherically symmetric (static) stars 
for models computed with stiffer EOSs. This is because such models are less centrally con- 
densed and gravitationally bound. This also explains why the momenta of inertia of rapidly 
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rotating neuron star configurations from the x = — 1 EOS show the greatest deviation from 
those of static models. 



4.3. Fractional moment of inertia of the neutron star crust 



As it was discussed extensively by lLattimer fc Prakashl (120001 ) (and others), the neutron 
star crust thickness might be measurable from observations of pulsar glitches, the occasional 
disrupts of the otherwise ext remely regular p ulsation from magnetized, rotating neutron 



stars. The canonical model of iLink et all (jl999l ) suggests that glitches are due to the angular 
momentum transfer from superfluid neutrons to normal matter in the neutron star crust, the 
region of the star containing nuclei and nucleons that have dripped out of nuclei. This region 
is bound b y the n eutron drip density at which nuclei merge into uniform nucleonic matter. 
Link et all (119991 ) concluded from the observations of the Vela pulsar that at least 1.4% 
of the total moment of inertia resides in the crust of the Vela pulsar. For slowly rotating 
neutron stars, ap plying several realist i c had ronic EOSs that permit maximum masses of 
at least ~ 1.6M© lLattimer fc Prakashl (2000) found that the fractional moment of inertia, 
AI/I, can be expressed approximately as 



AI 



28irP t R 3 (l- 1.67/3 -0.6/3 2 



3Mc 2 







1 + 



2P t (l + 5p- Up 2 

p t m b c 2 (3 2 



(22) 



In the above equation AI is the moment of inertia of the neutron star crust, / is the total 
moment of inertia, (3 = GM/Rc 2 is the compactness parameter, mj is the average nucleon 
mass, p t is the transition density at the crust-core boundary, and P t is the transition pressure. 
The determination of the transition density itself is a very complicated problem. Different 
approaches often give quite different results. Similar to determining the critical density for 
the spinodal decomposition for the liqu id-gas phase transit ion in nuclear matter, for uniform 
npe- matter, lLattimer fc Prakashl (120001 ) and more recently iKubis I (120071 ) have evaluated the 
crust transition density by investigating when the incompressibility of npe-matter becomes 
negative, i.e 



2 d E dE 2 
Kfj, — p + 2p—r- + 6 



dp 2 



dp 



P 



J2 r? 
2 u ^ sym 



dp 2 



+ 2p- 



dE 



sym 



dp 



- IE 



-l 

sym 



dE 



P- 



sym 



dp 



< (23) 



(see Fig. ^ where E (p) is the EOS of symmetric nuclear matter, E sym is the nuclear sym- 
metry energy, and 5 = (p n — p p )/( p v , + p p ) is the asymmetry parameter. Using this approach 
and the MDI interaction, IKubis I J2007I ) found the transition density of 0.119,0.092,0.095 
and 0.160/m -3 for the x parameter of 1,0,-1 and —2, respectively. Similarly, we have 
calculated the transition densities and pressures for the EOSs employed in this work. Our 
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Fig. 9. — (Color online) The incompressibility, K^, as a function of baryon density p. 
res ults are summ arized in Table HI We find good agreement between our results and those 



by iKubis I (120071 ) with the MDI interaction. It is interesting to notice that the transition 
densities predicted by all EOSs are in the same density range explored by heavy-ion reactions 
at intermediate energies. The MDI interaction with x = and x — — 1 constrained by the 
available data on isospin diffusion in heavy-ion reaction at intermediate energies thus limits 
the transition density rather tightly in the range of pt = [0.091 — 0.093] (/m~ 3 ). 

The fractional momenta of inertia AI/I of the neutron star crusts are shown in Fig. [TO] 
as computed through Eq. (1221) with the parameters listed in Table HJ It is seen that the 
condition AI/I > 0.014 extracted from studying the glitches of the Vela pulsar does put a 
strict lower limit on the radius for a given EOS. It also limits the maximum mass to be less 
than about 2M for all of the EOSs considered. Similar to the total momenta of inertia the 
ratio AI/I changes more sensitively with the radius as the EOS is varied. 



Table 4: Transition densities and pressures for the EOSs used in this paper. 



EOS 


MDI(x=0) 


MDI(x=-l) 


APR 


DBHF+Bonn B 


Pt(fm- 3 ) 


0.091 (0.095) 


0.093 (0.092) 


0.087 


0.100 


P t (MeV fm~ 3 ) 


0.645 


0.982 


0.513 


0.393 



The first row identifies the equation of state. The remaining rows exhibit the following quantities: 
transition density, transition pressure. The numbers in the parenthesis are the transition densities 
calculated bv lKubisI ((2003). 
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Fig. 10. — (Color online) The fractional moment of inertia of the neutron star crust as 
a function of the neutron star mass (left panel) and radius (right panel) estimated with 
Eq. (|22|) . The constraint from the glitches of the Vela pulsar is also shown. 

5. Summary 

Recent experiments in terrestrial nuclear laboratories have narrowed down significantly 
the range of the EOS of neutron-rich nuclear matter although there are still many remaining 
challenges and uncertainties. In particular, the EOS for symmetric nuclear matter was con- 
strained up to about five times the normal nuclear matter density by collective flow and par- 
ticle production in relativistic heavy-ion reactions. The density dependence of the symmetry 
energy was constrained at subsaturation densities by isospin diffusion and isoscaling in heavy- 
ion reactions at intermediate energies. Applying the EOSs constrained by the heavy-ion re- 
action data we have studied the neutron star momenta of inertia of both slowly and rapidly 
rotating models within well established formalisms. We found that the moment of inertia of 
PSR J0737-3039A is limited in the range of / = (1.30-1.63) x 10 45 (g cm 2 ). The neutron star 
crust-core transition density falls in a very narrow interval of pt = [0.091 — 0.093] (fm~ 3 ). 
The corresponding fractional momenta of inertia AI/I of the neutron star crust are also 
constrained. It is also found that the moment of inertia increases with rotational frequency 
at a rate strongly dependent upon the EOS used. 
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